function Fig4_e_f()
Fcolor='k';
Bcolor='w';
Fsize=18;  %font size

ED_pd_25 = readmatrix('processed data/df_ED_pd_25.csv');
nu = ED_pd_25(1, 2:end) - 1;
klb = ED_pd_25(2:end, 1);
map = ED_pd_25(2:end, 2:end);

f=figure();
s = pcolor(nu, klb, map);
s.EdgeColor = 'None';
ax = gca;
ax.YDir='normal';ax.XDir='normal';
f.Color = Bcolor;ax.Color=Bcolor;
ax.XColor=Fcolor;ax.YColor=Fcolor;
xl=xlabel('x');yl=ylabel('y');
yl.FontSize=Fsize;xl.FontSize=Fsize;
cMap=getColor('YlGnBu');
colormap(cMap);
ax.CLim=[-1.2,0.8];
axis([-0.85, -0.15, 0.5, 3]);


ED_pd_200 = readmatrix('processed data/df_ED_pd_200.csv');
nu = ED_pd_200(1, 2:end) - 1;
klb = ED_pd_200(2:end, 1);
map = ED_pd_200(2:end, 2:end);

f=figure();
s = pcolor(nu, klb, map);
s.EdgeColor = 'None';
ax = gca;
ax.YDir='normal';ax.XDir='normal';
f.Color = Bcolor;ax.Color=Bcolor;
ax.XColor=Fcolor;ax.YColor=Fcolor;
xl=xlabel('x');yl=ylabel('y');
yl.FontSize=Fsize;xl.FontSize=Fsize;
cMap=getColor('YlGnBu');
colormap(cMap);
ax.CLim=[-1.2,0.8];
axis([-0.85, -0.15, 0.5, 3]);



end